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ABSTRACT 

The potential of using interface layer to reduce thermal stresses in the matrix of 
composites with a mismatch in coefficients of thermal expansion (CTE) of fiber and matrix 
has been investigated. It was found that the performance of the layer can be defined by the 
product of the CTE and the thickness, and that a compensating layer with a sufficiently 
high CTE can reduce the thermal stresses in the matrix significantly. A practical procedure 
offering a window of candidate layer materials is proposed. 


*Visiting Professor from University Paris 6, LMT Cachan, 61 Avenue du President 
Wilson, 94230 Cachan, France 


INTRODUCTION 

Metal matrix composites reinforced with ceramic fibers are attractive because of their 
high specific stiffness and strength. Advantage can be taken from the high temperature 
Strength of ceramic fibers and the ductility of the metal matrix to produce a composite with 
superior combined properties. However, a thermal mismatch exists between the ceramics 
fibers and the metal matrix: the ceramics having a low coefficient of thermal expansion 
(CTE) and metals have higher values of CTE. This thermal mismatch induces stresses in 
the composite when subjected to temperature change. Matrix cracking has been observed 
after cooling down from processing temperature to room temperature for brittle matrix 
materials [1]. It has been proposed that the addition of an interface layer between the fiber 
and the matrix can reduce the tensile residual stresses in the matrix to a level which is low 
enough to avoid matrix cracking. 

Some numerical parametric studies have been performed [2], [3] in an attempt to 
determine the material parameters which are most beneficial in reducing thermal stresses. It 
has been suggested that the optimum interface layer should have a CTE between those of 
the matrix and fiber, with a low allowable modulus and a high layer thickness. A 
subsequent simplified analysis which retains the dominant physical features of the problem 
[4] showed that a compliant layer could not reduce the stresses in the matrix significantly. 
However, a layer with a high CTE, compensating layer, can reduce the stresses in the 
matrix significantly. The purpose of this study is to detremine the validity of the simplified 
method and to produce results which are useful for design. 

The study is based on a 3 cylinder model, isolating one fiber with an interface layer 
and a matrix layer (Fig. 1). Only monotonic cooling is studied and the variation of the 
materials properties with temperature is not considered. In a first part, the results of the 
simplified analysis which assumes a rigid fiber and a very thin layer, are recapitulated. 
This analysis identifies the properties of an interface layer. In the second part, the 
assumptions made previously are released and a complete elastic analysis is performed 
assuming that the three materials are isotropic and linearly elastic. A systematic parametric 
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study is also conducted and a practical procedure offering a window of candidate layer 
materials is proposed. 

In a later study [5] both the interface layer and the matrix are considered to be elastic- 
plastic, and the temperature dependence of the properties of the 3 materials is considered. 
In addition to the 3 cylinder model, a unit cell of a hexagonal array of the fibrous composite 
is defined and finite element computations are performed under several thermal-mechanical 
loadings. 
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NOTATIONS 


f 

fiber 

t 

interface layer 

m 

matrix 

j 

f, £, m 

R j 

external radius of j 

t/=R/-Rf layer thickness 

e j 

Young's modulus of j 

v j 

Poisson's ratio of j 

Xj and |ij 

Lamd coefficients of j 


EjV; E; 



V (l- 

2Vj)(l + vp ’ rj 2(1 + Vj ) 


CTE of j 

AT 

change in temperature 

a *j 

radial stress (<%) in j 

CT 9 j 

hoop stress (<799) in j 

^zj 

longitudinal stress (a zz ) in j 


radial strain (%) in j 

£0j 

hoop strain (£90) in j 

e zj 

longitudinal strain (e zz ) in j 


Mises equivalent stress in j 


The radii Rf and R m are related to the fiber volume 
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Cf by the relation 


1 . SIMPLIFIED ELASTIC ANALYSIS 
1.1 Stress-Strain Analysis. 

Ceramics fibers are usually four to five times stiffer than the metal matrix materials. 
Hence, in order to simplify the calculations, it is assumed that 


» 1 and 



The fiber volume fraction, Cf is assumed to be of the same order as the matrix volume 
fraction so that 


-^=1 

1-Cf 

For the given assumptions the stiff fiber controls the thermal expansion in the axial 
direction of the cylinder model. Hence, the stress distribution in the matrix and interface 
layer is governed by a plane strain problem with a fixed inner radius at the fiber interface. 
The analysis reveals that the differential CTE's of matrix and fiber, and interface layer and 
fiber dictate the thermal stresses. Consequently, the differential CTEs are defined as 

Aa m = a m - af and Aa i = 07 - ctf 


The stress distribution in the matrix is given by the plane strain solution for a thick 
walled cylinder subjected to an unknown internal pressure p, a traction free outer surface, 
and a temperature change AT. Superimposing the solutions given in [6] for an internal 
pressure and a temperature change gives the stress distributions in the matrix 


^rm — P 


(Rn,/Rf) 2 -' 
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The highest stresses in the matrix occur at the inner surface of the cylinder, r - Rf. The 
von Mises equivalent stress 


o = ^]^+gI + gI-g i g b -g i g z -g q g z 


can be calculated as 
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Cf 


1-Cf 


Cf 


[3 / Cf + 1 - 4v m (l - v m )] - pH m Aa m AT— —2[1 - 2v m ] 
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The layer is subjected to a pressure in the radial direction resulting from the contact with the 
matrix 


g t ( = -p 


The temperature change AT introduces stresses in the hoop and longitudinal directions due 
to the constraint from the fibers. The stresses in the hoop and longitudinal directions are 
equal and are the sum of the stresses caused by the pressure and the thermal expansion 


g Q( = g z t~ ~T~T 

\-v ( 


P + 


— - — E/Aa^AT 
1-v, 


6 


The von Mises equivalent stress in the layer is given by 


Of = — - — |E;Aa^AT + (1 - 2v/)p| 

1- v / 

Compatibility in the radial displacement at the boundary between the inner surface of 
the matrix and the interface layer is the condition which isolates the value of p which is 
found to be 


[l + v m ]- 


P = 


tj^Of l + v/ 
RfActml-V/ 


v 2 \ 
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1.2 Sensitivity Stu dy 

1.2.1 Stresses in the matrix 

The term E m Aa m AT enters in all the expressions for stress and all results are 
normalized with respect to this term. The stresses at the inner radius of the matrix are 
shown in Fig. 2 as a function of normalized pressure for a fiber volume fraction Cf = 0.4 
in a Ti 3 Al matrix with v m = 0.25. It can be seen that a pressure at the interface causes 
increased tensile stresses in the hoop and axial direction that could cause matrix cracking, 
initiated from defects at the fiber matrix interface. 

The value of the critical stress for matrix cracking, a cm , sets an upper limit for the 
principal stresses Ge m , a zm , and a rm and defines bounds on the pressure p if matrix 
failure is to be avoided. The allowable equivalent stress, Gy m , in the matrix defines other 
bounds on the pressure. These conditions define a window for the interface pressure, Fig. 
2, so that the failure criteria are not violated. 

Inspection of the expression of p reveals that the pressure may be reduced by the 
introduction of a compliant layer with low modulus which would increase the denominator, 
but the second term in the denominator of p is usually sufficiently small to be neglected in 
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Aa/t/ 

the calculations and the pressure is then linearly dependent on the parameter . This 

Aa m R f 

E t/ 

effect is demonstrated in Fig. 3 where it can be seen that the influence of m - is small. 

E * R f 

Figure 2 shows that the Mises equivalent stress has a minimum for an interface pressure 
that is close to zero. The pressure p is zero when 

A<x m R t m 1 + v, 

from which it can be deduced that the stresses in the matrix are strongly influenced by the 
parameter Aa^t/. The results from Fig. 2 and Fig. 3 can then be combined in Fig. 4 
where the highest principal stress and equivalent stress at the inner radius of the matrix 
cylinder are plotted versus the parameter Aa/tj/Aa m Rf. 

If the pressure is to be reduced to a minimum, the differential CTE for the layer has to 
be Aa/ ~6Aa m , when t//Rf = 0.1 For the case of SiC fibers, otf ~ 5 x 10‘ 6 /‘C in a Ti3 A1 
matrix, a m ~ 10 x 10' 6 , requires that a* ~ 35 x 10‘ 6 . Handbook values [7] of metals with 
a reasonably high melting point and high CTE are: silver 26 x 10 6 , hafnium 500 x 10' 6 
and copper 17 x lO' 6 . It is evident therefore that an interface layer of readily available 
materials with high CTE has the potential of substantially reducing thermal stresses in the 
matrix. 

1.2.2 Stresses in the Interface Laver 

The tradeoff for the improved stress state in the matrix is that the thermal mismatch 
has to be taken up by the interface layer. The deformation of the layer is constrained in the 
hoop and longitudinal direction so that tensile stresses develop in the layer. The highest 
principal stress in the interface layer is shown in Fig. 5 as a function of the different 
dimensionless parameters defining the problem. The stress is strongly dependent on the 
layer modulus, E^/E m , and CTE, Aa^/Aotm, but weakly dependent on the layer thickness 
t//Rf. This indicates that the stress is dominated by the constrained thermal expansion, 
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which is given by the second term in the stress expression, and is weakly dependent on the 
interface pressure, the first term in this expression. The product Aa/t/ governs the 
reduction of stress in the matrix, Fig. 4, and has to have a required value to reducing the 
stresses in the matrix to an acceptable level. It can be deducted from Fig. 4 and Fig. 5 that 
it is more favorable to have a thick layer and a moderate high layer CTE than a thin layer 
and a high layer CTE in order to fulfill the requirement of stress reduction in the matrix and 
to avoid high tensile stresses in the layer. 

2. ELASTIC ANALYSIS 
2.1 Stress-strain Analysis. 

The composite is not subjected to transverse loading and the outer surface of the 
compound cylinder is traction free. The thermal loading is axisymmetric with respect to the 
z-axis, then the displacement in the transverse plane is radial: Uj(r). The fibers are long 
and the stress and strain distributions are constant in the z-direction except at the end 
regions which will not be studied here. A generalized plane strain assumption is made: E z j 
= constant = e z . The other strains are given by: 

Uj 

e r j = Uj and e e j = — *■ 


and there are no shear strains. 

The linear thermo-elastic behavior gives the stresses as 


a n = + e z j + 2|ijUj - (3Xj + 2p.j|a jAT 

CT 0j = \/uj + ^ + e z ) + 2pj^ - (3Xj + 2pj)ctjAT 
r / r 


CT zj — 


UjH — ^ + e z + 2jije z — AT 

r 
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The equilibrium equations in cylindrical coordinates give: 


+ “( a rj - a 9j) “ o 

which gives the displacement field as 



Uj = Cjr + Dj— ■*- 


where cj and Dj are constants to be found by the boundary conditions. We have Df = 0 in 
order to obtain a finite solution in (f). 

To summarize, the strains are given by 


( Ri' 


e r j - c j 


D j 


fRA 2 


e 0j = Cj + 


V r J 


D; 


e z j = e z = constant 


The stresses are defined by 


fRA 2 

c r j = 2(Xj + |ij)Cj + A.je z -2p.j — Dj-(3X.j + 2^j)ajAT 


V r J 


fRA 2 

cjgj — 2(A.j + jij)cj + A.je z + 2jij — Dj - (3A.j + 2\i j)aj AT 

r ) 


a z j = 2XjCj + (>.j + 2pj)e z -(3Xj + 2|ij)ajAT 
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The strains and the stresses are uniform in the fiber. The axial strain is constant in the 
whole composite and the axial stresses are uniform in (£) and (m). The Mises equivalent 
stress Cj is defined by 

®1 (°rj- °8j) 2 + (o 8 j “ °zj) 2 + (°rj - °zj ) 

= 4i3f5i) 4 D? + ( Cj -e z ) 2 

_ V r / 

The six unknown Cf, C(, D(, c m , D m and e z are determined by the following conditions: 

Uf = U/ at r = Rf 
U/ =U m atr = R / 
a rf = a r/ at r = Rf 
c rl ~ a rm at r = R/ 

°rm = 0 at r = R m 

2tt | a z jrdr = 0=> Rfa z f + (R^-Rfja z /+(Rm-R/)ct zm = 0 
0 

The linear system obtained is given in Appendix I. 

Since ctrm(Rm) = the stresses in the matrix can also be written as 

CT 0m = 1 + ^ r j 
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CT zm ~ c m + e z ] 


It is reasonable to assume that D m > 0 (because this gives a compressive 0 ^ and a tensile 
O 0 m ). In this case, it is obvious that the minimal hoop stress is obtained for D m = 0, 
which gives = 0 and C 0 m = 0, i.e. a uniaxial long itudinal stress state. This result was 

confirmed numerically (see Section 2.3). 

2.2 Sensitivity Study 

The properties of the interface layer are defined by the four parameters: E (, V(, a.£ 
and i(. A sensitivity study was conducted in order to measure the influence of each of 
these parameters on the reduction of the residual stresses in the matrix. As a conclusion to 
this study, a procedure allowing the choice of candidate interface layer materials is 
proposed in Section 2.3. 

The fiber and the matrix are defined by the following data 


• Fiber: SiC (SCS 6 ) 

Ef = 360 GPa, Vf = 0.17, ctf = 4.9 x 10' 6 /*C 


• Matrix: TC 3 AI 

E m = 75.2 GPa, v m = 0.25, a m = 11.7 x lO'fy'C 


The radii Rf and R m are related to the fiber volume fraction by 


Cf = 



V^m ) 


= 40.5% 


The characteristics of the fiber and the matrix were kept constant but the parametric study 
can be readily applied to fibers and matrix with other properties. 

The range of interface layer parameters were chosen as follows: 


12 


69 < E^ £ 517 GPa 
4 < cc^ £ 30 x 10' 6 /*C 

0.028^ — <0.28 
Rf 

The Poisson's ratio of the layer vg was kept constant: vg = 0.3. 

To study the impact of the interface layer, a reference state was defined which 
corresponds to the state of stresses in the matrix in a 2 cyclinder without the layer. The 
stresses in the 3 cylinder model have been normalized with the reference stresses. The 
stresses in the layer and matrix are given at the inner radii. 

2.2. 1 Influence of the Young's Modulus Eg {ag, ig fixed) 

The results of the sensitivity study have shown that in the range 69 + 200 GPA, the 
value of E^ has little influence on matrix stresses (Fig. 6a). If Eg = 200 + 517 GPA, then 

if ag < 10 x 10'fy’C, the value of a m is again independent of the value of E^ 

(Fig. 6a). 

if ag > 10 x 10'fy‘C, the increase of Eg may reduce the stresses in (m) but 
increase the stresses in (£ ) to unacceptable values (see Fig. 6a and 6b). 

For the remaining of this study, we will choose Eg < 200 GPa. This is not too 
restrictive because an important property of the interface layer is a high value of CTE (see 
discussion in section 1.2) so that the layer is likely to be a metal for which the Young 
modulus is usually lower than 200 GPa. 

2.2.2 Influence of the CTE a g (Eg, t g fixed) 

The results show that generally the matrix equivalent stress decreases 
considerably as increases. However as illustrated in Fig. 7, the stresses reach a 
minimum after which they increase again with increasing ctg. 

The decrease of a m is more pronounced if the fixed thickness t g is high as is evident 
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from Fig. 7. 

In all the cases -^ DL < 1 : the addition of the interface layer does reduce the residual 
^ref 

stresses in the matrix. 

2.2.3 Influence of the Thickness t/ (E/, a/ fixed) 

The results are qualitatively the same as for CLf, i.e. in general a m decreases as \£ 
increases (see Fig. 8). 

2.2.4 A Convenient Parametric Representation 

The conclusions of the above sensitivity analysis are in agreement with those of the 
simplified analysis [4] described in section 1. That analysis lead naturally to the graphical 
representation given in Fig. 4. In the analysis it was shown that the important layer 
parameter which affected the matrix stress is Aa/t/. This graphical representation is used 
to present the results of the sensitivity study and they are given in Fig. 9. These results 
were generated using E^ = Em = 75.2 GPa, and they are nearly the same for all values of 
E£ lower than 200 GPa. The representation proves to be very useful since all the factors 
influencing the matrix stresses can be presented in one graph. Also shown in Fig. 9 are the 
results of the simplified analysis. The complete analysis predicts lower stresses than those 
from the simplified analysis, which is not surprising since the full analysis avoids 
assumptions which imply constraint. 

The range of x / = was limited to the domain of allowable pressure (defined 

Aa m R f 

in Fig. 2 and 3). The minimum value of a m corresponds to \£ ~ 0.65. If > 0.65, it 
was found that the radial stress in (m) becomes tensile and the hoop stress compressive 
(arm > °. <*8m < 0, a zm > 0). Layers such that x* > 0.65 are over-compensating layers. 

2.3 A Convenient Procedure & Choose Candidate Laver Materials 

The interface layer is defined by 4 parameters: E/, v^, a( and t^. In the sensitivity 
study we kept the Poisson's ratio of the layer constant: v/ = 0.3 and we found that it is 
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reasonable to take the Young modulus of the layer in this range: = 69 + 200 GPa. So 

we kept E/ = 69 GPa = constant. 

The matrix stresses can be readily calculated using Fig. 9. However, using this form 
requires further effort when making selection of the interface layer. This procedure is 
developed to deal with this situation. The axes used in the proposed graphs are the 
thickness of the layer and the CTE of the layer. The ratio of the matrix Mises equivalent 
stress with and without the layer is determined, and the graph shown in Fig. 10 then 
presents contours of constant reduction factor. Similarly in Figs. 1 1, 12, 13 the reduction 
factors for the damage stress, tensile hoop stress and tensile axial stress are plotted in an 
identical manner. 


Reduction of the Mises equivalent stress: 

• Of O / 

It has been noticed that ■=- L - and rr 1 - decrease if a.( decreases. So, for a given 

CT ref ^ref 

stress reduction inthe matrix (i.e. a given contour in Fig. 10) one can choose such that 

does not exceed a certain value to avoid layer cracking. The min. value of was 
^rcf ^ref 

found to be ~ 0.387 (see shaded line in Fig. 10). This contour was also found to 
correspond to a uniaxial longitudinal stress state in the matrix (see remark at the end of 
section 2.1). 


Reduction of t he Damage Equivalent Stress 
In the same way that the Mises equivalent stress is related to the distorsion energy, a 
so-called damage equivalent stress [8] is related to the total elastic energy and is defined by 


a* = G 


-(l + v) + 3(l 


-2v)f^f 

U; 


il/2 


where CT H = ^a kk . 
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One notices the presence of the triaxiality ratio which is known to be a main feature in 

(T 

Failure. 

The contours of Fig. 11 can be used in the same way as those of the Mises stress. 

* ♦ 

The same remark about the decrease of Of and O/ when (X/ decreases applies also. If we 
compare the contours of the Mises and the damage equivalent stresses, we find that: 


i) In the "upper" region of the plane defined by a/ = 30 x 10 -6 /°C, — = 0.28 

Rf 

* 

and — iP- = 0.387 (shaded line) the criterion based on the reduction of the Mises stress is 
°ref 


f 


a bit more severe than the second one, i.e. for a given point 
q* o 

CT ref a ref 


1L 

Rf 




we have 


ii) In the remainder region, the criterion based on the reduction of the damage 

equivalent stress is a bit more severe than the first one. The reason is the following. 

* 

The calculations have shown that a re f = a re f , so 


£ H <1 

a 3 

so, the triaxiality ratio in the matrix is 
£ j in the "upper" region i) 

> — in the remainder region ii) 



Reduction of the Hoop Stress 

It appears that the reduction of the hoop stress in (m) is much more important than the 
reduction of the Mises stress. The hoop stress can actually be decreased to a zero value, 
which is shown by the shaded line in Fig. 12. This line corresponds also to the min. Mises 
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stress and to a uniaxial z-stress in the matrix. In the "upper" region defined by the shaded 
line as previously, the hoop stress in (m) becomes compressive: O0 m < 0. This region 
corresponds also to x/ > 0.65, as defined in section 2.2.4. 

Reduction of the Longitudinal Stress 

It appears that the reduction of is smaller than the reduction of a m or <70 m . The 

shaded line (o m minimal and G0 m = 0) gives = 0.6. A lower reduction cannot be 

a zref 

expected and if the reduction of a zm is a concern, then the window of candidate layer 
materials is reduced as shown by comparing Fig. 13 to the previous figures. 

CONCLUSIONS 

The thermal mismatch between the fiber and the matrix has to be taken up by the 
interface layer and can subject it to high stresses if the Young modulus of the layer is high. 
If this modulus is taken in a certain range, common to most metals, then the layer 
performance is defined by the product of the CTE and the thickness. 

A compensating layer with a sufficiently high CTE has the potential of reducing the 
thermal stresses in the matrix significantly. Both the CTE and the thickness of the layer can 
be adjusted in order to keep the stresses in the layer under a certain level. The maximum 
hoop stress in the matrix can be reduced substantially but the axial stress in the matrix is 
less affected by a layer. This implies that compensating layers can be expected to be very 
successful in preventing cracking in composites where predominantly radial cracking is 
observed in the matrix, but that while axial stresses can also be reduced, this reduction is 
less dramatic. 

Graphs have been produced which are simple to use. From Figs. 10 to 13, the 
reduction in stress for any value of layer CTE and thickness are readily determined. The 
absolute values of matrix stress can be easily determined by reference from another diagram 
(Fig. 9). 
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APPENDIX: Linear System for the Constants of the Elastic Solution 
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Fig. 2 Stresses at the inner surface of the matrix cylinder as a function of normalized 
pressure. 
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Fig. 3 Normalized pressure as a function of the layer CTE and thickness. 
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Fig. 4. Maximum principal stress and Mises equivalent stress at the inner surface of the 
matrix cylinder as a function of the layer CTE and thickness. 
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Fig. 5 Longitudinal and hoop tensile stress in interface layer. 
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Fig. 6(a) Influence of the layer Young modulus on the stresses in (a) matrix, and (b) 
layer. 
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Fig. 6(b) Influence of the layer Young modulus on the stresses in (a) matrix, and (b) 
layer. 
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8 Influence of the layer thickness on the stresses in the matrix. 
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